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ABSTRACT 



The Hydro-Numerical Prediction Model developed by Hansen 
is applied to San Diego Bay, and the results compared both 
with the hydraulic model and the real data obtained by 
field measurements. This allows one of the few good compar- 
isons between •■numerical and hydraulic models for the 
prediction of actual conditions. 

The bay was divided into two sections that were run 
separately in order to obtain the desirable spatial resolu- 
tion. This division required solving the problems of proper 
tuning and matching techniques between both portions. The 
solution involved the addition of an appended pseudo-bay 
to the first section of the model in order to compensate 
for the correct tidal prism. The effects of a proposed 
second open entrance in the southern part of the bay were 
studied. This resulted in an increase of flushing in the 
southern portion of the bay but caused the currents in the 
center of the bay to be small which decreased dispersion 
in the central portion of the bay. In general, both models 
produced similar and reliable results, but there was a 
considerable reduction of cost and time with the numerical 
model. 
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I. INTRODUCTION 



A . REVIEW 

The study of bays, estuaries and other semi-restricted 
areas is difficult because these masses of water are subject 
to constant movement with short time irregular periodicities. 
The important mechanism driving the circulation in these 
areas generally is the tide. The water is partially re- 
newed and mixed during each tidal cycle due to the inter- 
change of water between such areas and the ocean. The 
shape and depth of the area governs the direction and speed 
of the currents generated during this interchange. The 
constant flow of water causes a constant state of non- 
equilibrium of constituents, temperature and potential 
energy inside the bays or estuaries. 

The study of such areas by field observations is ex- 
tremely laborious requiring large expenditures of manpower 
and money. Another approach to the solution of the problem 
is to reproduce the physical characteristics of these areas 
by the use of hydraulic models. Hydraulic models have been 
used for many decades and have proven to be a reliable 
method for studying an area although space and operating 
problems make them a relatively expensive tool of 
investigation. 

With the advent of high speed computers, the numerical 
solution of the differential equations used to describe 
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flow was made possible. Oceans, seas, gulfs and other 
water masses have been described by the use of this tech- 
nique. It transforms a physical model into a program 
which can be modified to simulate different conditions more 
easily, can be stored for future uses and alterations, is 
faster and allows solutions to be obtained at reasonable 
...cost . 

The U. S. Army Corps of Engineers made a study of San 
Diego Bay by means of a hydraulic model based on field 
measurements made in early 1967. The results have been 
reported in a Technical Report [10] and published in a 
paper [1] by the U. S. Corps of Engineers. 

The Hydro-Numerical Prediction Model developed by 
Hansen is applied to San Diego Bay, and the results com- 
pared both with the hydraulic model and the data obtained 
by field measurements. This allows one of the few good 
comparisons between numerical and hydraulic models for 
the prediction of actual conditions. 

B. OBJECTIVES 

The objective of this thesis is to reproduce by means 
of the Hansen hydrodynamic -numerical model the existing 
conditions in San Diego Bay and compare them with condi- 
tions produced in the hydraulic model operated by the U. S. 
Corps of Engineers. This comparison will result in a 
proper evaluation of the merits and limitations of both 
models . 
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The division of the bay into two separate models is 
attempted in order to obtain comparable resolution in the 
numerical model. This division requires solving problems 
in the method of running numerical models in series and of 
transfering boundary conditions from one section to 
another. 
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II. HYDRAULIC MODEL OF SAN DIEGO BAY 



A. DESCRIPTION OF SAN DIEGO BAY 

San Diego Bay is located on the southern coast of 
California, U.S.A. The bay is long, curved, has an area 
of approximately 90 sq. km. and a maximum width of 3-1 km. 
[Fig. 1]. It is connected with the ocean .by the narrow 
Zuniga Channel at its northern end. 

The inflow of fresh water can be considered negligible 
except during local heavy rainfall. The circulation of the 
entire bay is driven by the tide at its unique entrance. 

The tidal circulation can be perturbed by local wind 
conditions . 

The tide is of the mixed type with a period of 24.83 hours 
and maximum range of about 1.89 meters. Field measurements 
indicate that tidal variations become larger in the bay 
increasing to 1.98 meters at the center portion of the bay 
and 2.1 meters near the southern end. Currents have a 
maximum speed of about 0.7 meters/second in Zuiga Channel. 

A channel is maintained for navigation starting at 
Zuniga Channel and continuing throughout San Diego Bay. 

The channel depth is 12.8 meters below MLLW from the ocean 
to the vicinity of the South Bay gage. Other channels of 
minor importance exist in the southern portion of the bay. 
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FIGURE 1 
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B. PURPOSE OF THE U.S. CORPS OF ENGINEERS MODEL STUDY 

The natural shape of San Diego Bay and its only entrance 
at the northern end, its use as a Naval Base, and the in- 
fluence of a large community surrounding the Bay, have 
brought into consideration the construction of a second 
navigation entrance somewhere along Silver Strand Beach. 

This second entrance is expected to cause decreasing ship 
traffic throughout Zuniga Channel and increasing flushing 
rate of the Bay, possibly improving the environmental 
conditions within the bay. A major consideration in eval- 
uating the overall benefits of a second entrance depend 
heavily on assessing the extent to which it would increase 
the flushing rate. These effects could not be computed 
reliably by the available analytical methods [10], so a 
physical hydraulic model of the bay, in which several 
suggested locations for a second entrance could be tested, 
was constructed to investigate the effects of each proposed 
location, in detail. 

C. THE HYDRAULIC MODEL 

The dominant forces in rivers, narrow estuaries and 
some lakes, are pressure gradient forces driving circula- 
tions which are opposed by friction [12] . In large scales 
of circulations, the opposing forces are combinations of 
friction and Coriolis force, or Coriolis force alone. 

Similarity can be considered as being graded from 
geometrical through kinematic to dynamic in the attempt to 
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reproduce the existing conditions in a scale model. A 
model has geometric similai’ity if the ratios of all homol- 
ogous dimensions are equal. The model is said to have kine- 
matic similarity if the paths and patterns of motion are 
geometrically similar to those of homologous occurrences 
in the prototype, and the ratios of homologous velocities 
are equal at all times. Dynamic similarity is achieved 
when the ratios of homologous masses and forces affecting 
motion are equal at all times. Complete dynamical similar- 
ity is difficult to achieve in fluid models and, kinematic 
similarity, of necessity, involves elements of dynamic 
similarity . 

The general equation of motion can be written in the 
form: 



p^r = -p2Px^ - Vp + plcg + P 



Where : 




inertial force 



p2^x^ = coriolis force 



VP 



pressure gradient force 




gravity force 



F 



friction force 



Each term of the equations of motion represents the 
forces present in a prototype and what is to be represented 



in a model. Since the totality of forces cannot be easily 
reproduced in hydrodynamic models, a few important terms 
are usually singled out. 

In scaling a model, each quantity and dimension of those 
particular terms of the equation of motion used to describe 
a flow must be reduced to the model scale. Dimensionless 
ratios of terms are used to achieve correct model scales. 

The numerical value of the dimensional ratio must be the 
same for the model as for the prototype for simiarity to 
prevail . 

Several dimensionless ratios of terms are common to 
oceanographic modeling, and each one is used depending on 
the characteristics of the flow in the prototype. The 
Proude number is generally used in models of estuaries 
because it involve the ratio of inertial to gravity forces. 
These forces are, with friction, the dominant forces. The 
direction of flow is primarily down the gradient of 
pressure . 

The Froude number is defined as : 

p r _ inertial term _ V c 

gravitational term , /0 

(Ly/p) 17 ^ 

Where : 

V c = characteristic velocity 
L = characteristic length 
y = specific weight 
p = fluid density 
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As a preliminary step in scaling a phenomenon down to 
model size, it is useful to construct a ratio of units which 
by convention is: 




where is a number representing the dimensions of a cer- 
tain property of the model and A the dimensions of the 
homologous property in the prototype. 

The length ratio is fixed first in order to match the 
available space or equipment. The numerical value of each 
dimensionless group constant can be held altering the ratios 
of the other units included in a dimensionless group. By 
these procedures, time and other unit dimensions of other 
related properties may have to be altered to permit the 
fluid motion in the model to be regarded as behaving in 
the same physical manner as that in the prototype. 

The ratio of inertial (velocity and length) and gravity 

is of concern for the Proude Number. Shallow-water gravity 

waves, such as tides, propagate at a speed C = ^ gh. 

The velocity ratio (R ) is satisfied using the length 

c 

ratio (R t ) , where 
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the length ratio and velocity ratio are left as adjustable 
parameters 



It can be seen from this equation that R /( _ shortens with 

increasing and lengthens with increasing R^. 

A distinction between the vertical length (h) and the 

horizontal length (L) is often made. If the Proude model 

is geometrically similar to its prototype, R^ = R . 

The vertical scale (R^) generally must be chosen so 

that the least depth of water in which the flow is to be 

studied is approximately one cm. to prevent capillary 

effects from being important. Then, it may be necessary 

to depart from similarity by distorting vertical dimensions 

such that R, > R T . One is led to a choice of scales by a 
h L 

path which touches first the minimum depth requirement and 
then an adjustment of the R T and R such that either one 

Li U 

becomes constant, leaving the others to be adjusted. 

The R^ (time ratio) is usually the most flexible 
parameter available and determines the frequencies of the 
simulated tides. 

The effects of wind stress on the circulation in estu- 
aries and bays are often important and sometimes dominate 
the tide. To scale it in a model requires more knowledge 
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of the dynamics of momentum transfer across the air-sea 
interface than is now available. Empirical procedures 
must be made to suffice. 



For the U.S. Corps of 

following characteristics 

Area: 280 sq. km. 

Horizontal scale (R^ ) 

Vertical scale (R, ) 

h 

Velocity scale (R ) 
Time scale (R^) 
Discharge scale (R 
Volume scale (R ) 
Physical dimensions 



Engineers hydraulic model, 
and scales were used: 

= 1:500 
= 1:100 
= 1:10 
= 1:50 
= 1:500 000 
= 1:25 000 000 
= 35 x ^0 meters 



the 



D. RESULTS OF THE HYDRAULIC MODEL 

The type of flow being studied determines the predomin- 
ant forces to consider for dynamic similarity. The San 
Diego Bay model used a Froude number scaling criteria. 

After determining the dynamic scaling, a hydraulic model 
is constructed to match the scaled geometry of the prototype. 
After reproducing the geometry as closely as possible, kine- 
matic similarity must be obtained. This is accomplished 
by matching the hydraulic head and flow patterns at various 
locations in the model over the tidal cycle. This matching 
is generally accomplished by using roughness elements which 
can consist of concrete blocks or metal strips. The rough- 
ness elements have the effect of increasing turbulence and 
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changing mean flow patterns. This is a rather laborious 
and tedious procedure and usually is the most time consuming 
aspect of hydraulic model building. After the kinematic 
similarity is satisfactorily accomplished, the model is 
presumed to be calibrated. 

Tests of existing conditions were made under carefully 
controlled conditions of tides, currents, and simulated 
pollution input [1]. The results of these calibration 
tests, called Base Tests, were used to evaluate the effects 
of navigation openings. Thus, any differences noted during 
the tests of the proposed second entrances were attributed 
to influences of the plan being tested and not to errors 
in the construction of the model. 

The wind conditions on the day of observations can in- 
troduce large modifications to the normal surface and cur- 
rent characteristics in a bay such as San Diego Bay. Dif- 
ferences from day to day in the values of a mixed tide make 
a one-day observation a non-periodic function. Therefore, 
modifications of the field measurements were introduced in 
order to make them periodic and comparable to a non-wind 
situation . 

The effects of the south entrance on tides throughout 
the bay are shown by comparative tide curves for measured, 
base and planned conditions in figure 2-1. The second 
entrance caused a reduction in the tide range at the southern 
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HEIGHTS IN THE PROTOTYPE AND IN 
THE HYDRAULIC MODEL 



FIGURE 2-1 



23 



end of the bay by lowering high water and raising low water. 
At the south bay gage, the reduction in tide range was about 
15 cm. 

The effects of the second entrance on mid-depth current 
velocity are shown in figures 2 and 3- The opening reduced 
the maximum mid-depth current velocities at all stations on 
gauges "A" through "D" by 15 to 75 cm. /sec. The flushing of 
the southern portion of the bay became governed by the new 
entrance and there seemed to be little interchange of water 
between the northern and southern portion of the bay with 
the area in the vicinity of Pier No. 2 gauge having very 
small currents. Under those new conditions, the entire bay 
appears to have two completely different circulation 
systems driven by tides at Zuniga Channel and at the second 
entrance. The flushing appears to be increased in the 
southern portion by including a second entrance. 

The pictorial representations of the flow in the bay 
under both conditions [1] are shown in figures 29, 30, 36 
and 37. In these figures, the flow is represented at 
different stages of the tide, and its intensity is scaled 
by the length of the lines. 
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III. THE HANSEN HYDRO DYNAMICAL MODEL 



The assumptions of the numerical model are: 

1. The fluid is homogeneous and incompressible. 

2. Pressure is hydrostatic and thus, the changes in 
pressure are due solely to changes in water 
surface elevation. 

3. Advection terms are ignored. 

*1. The fluid is in hydrostatic equilibrium in the 
vertical direction. 

5. The geographical and vertical variations of the 
Coriolis force are neglected. 

Applying these assumptions to the equations of conserva- 
tion of momentum and mass, basic equations are obtained for 
the single-layer model developed by Hansen [6]. Horizontal 
momentum equations are integrated over depth to give 




6 fx + Kv2u “ xb ( x ) + T ( x ) 




g 1^- + Kv 2 v - T b (y) + t (y ) 




0 



The wind stress components are represented by 



t (x) 





27 



and the bottom stress components are represented by 



The various terms of the equations are defined as: 



K = coefficient of horizontal kinematic eddy 
viscosity 

X = drag coefficient 



n = surface elevation 
H = total depth = h + n 
x,y = space coordinates 
u,v = velocity components 
g = acceleration of gravity 
f = coriolis parameter 



Both the wind stress and bottom stress equations are 
empirically developed. The wind stress terms, as used in 
the model, are assumed to be a quadratic expression in wind 
speed, where X is the wind drag coefficient with a typical 
value of 0.65. The bottom stress is assumed to be non- 
linearly dependent on u and v, and its formulation, like 
the wind stress, have been derived by empirical means. It 
was originally formulated for shallow water application. 
Over deep water, its value becomes small and its 
applicability is questionable. 



b, x 

t (x) 




and 




j = wind speed component 
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An implicit central difference scheme is used for 
achieving time dependent solutions of the equations. At 
the boundaries, the values of n, u and v are taken at the 
actual points rather than from the surrounding points as 
shown in Figure *1 . 

The driving forces are tides input at the open boun- 
daries, and wind at the surface over the entire grid. The 
tide values are computed at each time step and introduced 
as new values of n at each point of the boundaries. 

The finite approximations are given below. The water 
surface elevations are first calculated using the 
conservation of mass equation. 

t+— t — — 

n 2 (n,m) = n 2 (n,m) - (H^ (n ,m)U fc (n,m) 

- H (n,m-l)U t (n,m-l) + (n-1 ,m) V t (n-1 ,m) 

- H v fc (n,m)V t (n,m) } 

The horizontal and vertical velocity components are then 
determined from respective horizontal and vertical momentum 
equations . 
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U t+At (n,m) = {1 - [At r/Hu t + At (n,m) ]\/u fc (n ,m) 2 +V (n,m) }U fc (n,m) 

t+ At 

+ At f V ^(njin) - n 2 (n,m+l) 

r) t+ ~(n ) m)} + At X t+At (n,m) 

t-i--- / 

V t+At (n,m) = {1 - [At r/Hv 2 (n ,m) ]Jv^ (n,m) 2 +U ^ (n ,m)} (n ,m) 

t+— 

- At f U ^ (n ,m) - - {j - ' 1 ( n 2 (n,m) 

t+— 

- n 2 (n+ljin)} + At Y^^Cnjin) 

where 

U^(n,m) = aU^Cn^m) + — (n-1 jirO+U^ (n+1 ,m) +U^ (n ,m+l ) 

+U fc (n,m-1) } , 

V^(n,m) = aV^(n 3 m) + ^— { (n-1 3 m) +V^ (n+1 ,m)+V^ (n ,m+l ) 

+V fc ( n , m- 1 ) } , 

U k(n,m) = (n ,m-l )+U t (n+1 jin-1 )+U^ (n 3 m)+U^ (n+1 jin) } , and 

V ^(n 3 m) = { V 1 " (n , 111 - 1 ) +V^ (n+1 > m-l)+V^' (n,m)+V^ (n+1 ,m) } . 



The factor a can be interpreted as related to the horizontal 
kinematic viscosity parameter. In the program it is treated 
as a tuning parameter and is related to the eddy viscosity 
coefficient by: 
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